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Abstract 



Double-diffusive convection driven by both thermal and compositional buoyancy in a rotating cylindrical annulus with conical caps 
is considered with the aim to establish whether a small fraction of compositional buoyancy added to the thermal buoyancy (or vice 
CN| versa) can significantly reduce the critical Rayleigh number and amplify convection in planetary cores. It is shown that the neutral 
Uh surface describing the onset of convection in the double-buoyancy case is essentially different from that of the well-studied purely 
Cj-ifhermal case, and does indeed allow the possibility of low-Rayleigh number convection. In particular, isolated islands of instability 
are formed by an additional "double-diffusive" eigenmode in certain regions of the parameter space. However, the amplitude of 
such low-Rayleigh number convection is relatively weak. At similar flow amplitudes purely compositional and double-diffusive 
£N| pases are characterized by a stronger time dependence compared to purely thermal cases, and by a prograde mean zonal flow near 
the inner cylindrical surface. Implications of the results for planetary core convection are briefly discussed. 
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(J planets is of significant interest as it drives the dynamo pro- 
" ^ cesses that ge nerate and sustain t he global ma g netic fields of 
i>~»ihese bodies ( Kono and Robertsl 20021 : JonesL 2007 ). Core 
convection is a double-diffusive process driven by density 
, variations due to non-uniform temperature and composition 
dBraginskv and Roberts . 1995 ). While ouble-diffusive phenom- 
ena are well-studied in oce anography, metallurgy, mantle con- 
vection and other contexts (IHuppert and Turnerl Il98lt iTurnei , 



1. Introduction 



Convection in the cores of the Earth and the terrestrial 



> 

o 



1974, 1985; Schmitt, 1994), their manifestations in core con- 



O 



vection remain poorly understood. It is thought that thermal and 
composi tional buoyancy in the Earth's core hav e comparable 
strength dLister and Buffettl[l995llNimmcJ.l2007l) . and that tem- 
perature and concentration of light elements have widely differ- 
ent molecular diffusive ti me scales, boundary conditio ns and 
i~ I source-sink distributions (IBraginskv and Robertsl 119951) . Yet, 
j>. most planetary and geo-dynamo models consider only thermal 
T convection or, at best, lump temperature and concentration into 
rN a single "codensity" variable. The last approach is poorly jus- 
g$ tified, as it is only valid for equal diffusivities and identical 
boundary conditions. Indeed, while eddy diffusivities due to 
small-scale turbulent mixing tend to attain similar values, the 
turbulence in many cases, e.g. weakly-convecting stratified lay- 
ers, may not be as fully developed to wipe out the la r ge dif - 



ferences in molecular diffusivities dBusse and Simitevl 1201 ll) . 
At the same time, relatively small variation s in diffusivity ratios 
may h ave significant dynamical effects (e.g. lSimitev and Busse , 



20051) . So far, only few studies have been published where 



thermal and composition al buoyancy are considered separately. 
Cardin and Olsonl ( 1992 ) performed an experimental investiga- 



tion of thermochemical convection in rotating spherical shell. 
A double-diffusive numerical dynamo model with a partly sta- 
ble thermal gradient and destab ilizing compositiona l compo- 
nent has been recently studied bv lManglik et al. (2010), as a sit- 
uation likely applicable to Mercury. Various driving scenarios 
where thermal and compositional grad ients are both destab iliz- 
ing have been explored numerically bv lBreuer et al.l (120101) . All 
of these papers report significant differences in their results to 
the single-diffusive (codensity) case and emphasize the need for 
further investigation. The onset of double-diffusive con vection 
in an ax isymmetric rotating system has been studied by Busse 
d2002bl) in certain asymptotic limits, and it was found that a 
small fraction of compositional buoyancy could significantly 
reduce the critical Rayleigh number, and thus amplify core con- 
vection. This prediction is potentially very important, as it may 
shed light on the thermodynamic state of the core and the en- 
ergy budget of the geodynamo. However, concerted numerica l 
simulations have so far failed to confirm it ( Breuer et al. , 20 lOh . 
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With this motivation, the goals of this letter are to establish 
the possibility of low-Rayleigh number double-diffusive con- 
vection, and to elucidate the mechanisms by which thermal and 
compositional buoyancy interact. To this end, a simple model 
of a rotating cylindrical annulus with conical end caps is con- 
sidered here. This model has been very useful in capturing the 
basic behaviour of nearly geostro phic convection in the equa- 
torial regions of planetary cores dBusse, l2002at IJonesi 120071) 
and offers significant mathematical and computational advan- 
tages. The attention is restricted here to the effects induced 
by the difference in diffusivity values, while the more realistic 
cases of distinct boundary conditions and source-sink distribu- 
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Figure 1: Sketch of the rotating cylindrical annulus with conical end caps. Note 
that the sketch is not to scale with the limits of a small gap and a small angle of 
inclination of the conical caps. 



tions are disregarded at present. The mathematical formulation 
and the methods of solution are presented in Section [2] Sec- 
tions [3] and [4] describe linear and finite-amplitude properties of 
double-diffusive convection. Conclusions and possible impli- 
cations for planetary cores are discussed in Section [5] 

2. Formulation and methods of solution 

A cylindrical annulus with conical caps full of a two- 
component fluid, and rotating about its axis of symmetry with 
an angular velocity f2 is considered. The configuration is shown 
in figure [T] an d a ma t hematical for mulation of the problem 
given earlier bv lBussd (119861 l2002bl) is adopted. In particular, 
the inner and outer cylindrical walls are kept at constant tem- 
peratures Tq + AT/2, and at constant values of the concentration 
of the light element Co + AC/2, respectively, such that a den- 
sity gradient opposite to the direction of the centrifugal force is 
established as the basic state of the system. The effect of the 
centrifugal force is similar to that of gravity in self-gravitating 
spheres and shells in that the buoyancy-driven motions occur in 
the same way as in the case when the force and the gradients of 
temperature and concentration are reversed. This formulation 
has the important a dvantage of being amenabl e to experimen- 
tal realizations (e.g. jBusse and Carrigan[|l974l) . The gap width 
d of the annulus is used as a length scale, d 2 /v - as the time 
scale, and vAT/k and vAC/k- as the scales of temperature and 
concentration of light material, respectively. Here v is the kine- 
matic viscosity, and k is the thermal diffusivity. A small-gap 
approximation, d/ro <K 1, is assumed, where ro is the mean ra- 
dius. This makes it possible to neglect the spatial variations of 
the centrifugal force, and of the temperature and concentration 
gradients of the static state, and to introduce a Cartesian system 
of coordinates with the x-, y-, and z-coordinates in the radial, 
azimuthal and axial directions, respectively. The Boussinesq 
approximation is adopted, in that the variation of density, 



p = po(l - y,AT(x - ®/F) - y s AC(x - T/L)), 



(1) 



is only taken into account in connection with the body forces 
acting on the fluid. Here, y, and y s are the coefficients of ther- 
mal and chemical expansion, and the other symbols are de- 
fined below. The linear dependence on x is unrealistic for the 
concentration as it requires zero-concentration boundary con- 
ditions ([5]). A more realisti c zero-flux condition would m ake 
the problem rather involved (IBraginskv and Robertsi Il995l) and 



will divert from the main focus of this paper which is to inves- 
tigate the influence of diffusivities isolated from the effects of 
boundary conditions. A discussion of different types of bound- 
ary conditions related to core convection and the geodynamo 
can be found, for exampl e, in ( Kutzner and Christensenl 2002 : 
Busse and Simitev, l2006h . Assuming a small angle of inclina- 
tion of the conical end caps with respect to the equatorial plane, 
and taking into account that the annulus is rotating, the velocity 
obeys approximately the Proudman-Taylor theorem and can be 
described in first approximation by its geostrophic part 



u = V x kt(/(x,y, t) + O(rjo), 



(2) 



where 770 «: 1 is the tangent of the said angle. Averaging over 
z, the governing equations for the leading order of the dimen- 
sionless deviations of the temperature ©, the concentration T, 
and the stream function x/i from th e static state o f no flow can be 
written in the 2D cartesian form (IBussel l2002bh 

{{d, - V 2 ) + J^)VV - rfdyip + d y (R t © + R s F) = 0, 

P(d, + J^)0- V 2 + d v tA = O, (3) 

p{d, + J»r - it 1 v 2 r + dyf = 0, 

where = (<3 v t/03 A - (d x iff)d y , and the definitions of the rota- 
tion rate, Prandtl, Lewis, thermal and compositional Rayleigh 
numbers 77*, P, L, R t , and R s are 
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respectively. Here, D is the diffusivity of the light material, h 
is the axial length of the annulus, and g = fi 2 ro is the average 
centrifugal acceleration analogous to gravitational acceleration. 
Fixed temperature and concentration, and stress-free BCs for 
the velocity are assumed at x = +1/2, 



^ = d 2 x i// = & = T = at 



1/2, 



(5) 



while periodicity is imposed in the y-direction. For further de- 
tails on the assumptions and for evidence of the utility of this 
model to capture the dynamics of convection in rotating spheri - 
cal sh e lls, th e reader is referred to the reviews of lBusse l (12002al) : 
Jonesl (120071) and the references cited therein. 

The linearized version of equations (|3} allows an analytical 
solution. The nonlinear equations (0 are solved numerically 
by a modi fication of the Galerkin spectral method used pre - 
viously bv lOrandBussd dl987h : ISch"naubelt and Bussd Jl992h . 
The dependent variables if/, and F are expanded in functions 
satisfying the boundary conditions 
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sinO(ji:+ 1/2)). 



After projecting equations Q onto the respective expansion 
functions, a system of nonlinear ordinary differential equations 
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Figure 2: The critical Rayleigh number i? ( cr jt of purely-thermal convection as 
a function of the radial wave number / for a = 5, if = 600, R s = 0, and 
P = 10~ 2 , 10~', 1, 10, 10 2 (from bottom to top). The solid circles indicate the 
preferred values of /. 



is obtained for the unknown coefficients ai„(i), ai„(t), bi„(t), 
bin(t), Q n (f) and ci„(t). The system is integrated in time by a 
combination of an Adams-Bashforth scheme for the nonlinear 
terms and a Crank-Nicolson scheme for the diffusion and the 
other linear terms. A truncation scheme must be introduced in 
practice: the equations and the corresponding coefficients are 
neglected when / > N x and n > N y , where the truncation pa- 
rameters N x and N v must be sufficiently large so that the physi- 
cal properties of the solution do not change significantly when 
their values are increased. The computations reported in the 
following have been done with /3 = 1, N x = 35 and N y = 55. 

3. The linear onset of double-diffusive convection 

Without loss of generality, small perturbations about the state 
of no motion can be assumed to take the form 



(if/, 0, F) T = (iff, 0, f ) T sin (ln(x + 1 12)) e 



iay+At 



(7) 



where a and / denote the azimuthal and the radial wave num- 
bers, A — cr + io), with cr e R and to e R being the growth 
rate and the frequency of oscillations, respectively. The su- 
perscript T denotes transposition, and (if/, 0, f ) T is a constant 
component vector. Then, the linearised version of equations (f3]) 
reduces to a matrix eigenvalue problem for A and (iff, 0, f ) T , 



a- 
a 
—I— 
P 

a 

~ l P 

;2„2 , 



. aR t . aR s 



a 




PL 



"d 2 


'V 














= A 





(8) 


a 2 











where a 2 = l l n z + a 2 . In the rest of this section attention is 
restricted to a single-roll convective structure in radial direction 
by setting 1=1. Equatorially-attached "multicellular" thermal 
convection has been previously found in s hells and annuli of 



finite gaps and convex ly curved caps (e.g. lArdes et al.L 11997 : 
Plaut and BusseL 120051) . However, these geometries are quite 



different from the small-gap limit considered here as they pro- 
vide radially inhomogeneous conditions for convection. Figure 
[^demonstrates that, in the small-gap limit and for the parame- 
ter values discussed below, / = 1 is always the preferred radial 
mode for the onset of thermal convection; the nonlinear results 
of section 2] further confirm the / = 1 assumption. 

The solution to problem ^ can be found in analytical form, 
and figure [3] shows the growth rate of the perturbations, cr = 
Re(A), as a function of the thermal Rayleigh number R, for 
fixed values of the other parameters. The eigenmodes of purely 
thermal convection are also shown in the figure for compari- 
son. Because the matrix in ^ is of size 3 x 3, it can have 
up to three distinct eigenmodes for typical parameter values. 
The analogous eigenvalue problem for purely thermal convec- 
tion has a matrix of size 2x2 that can have up to 2 eigenmodes 
at most. Thus, a basic distinction between purely thermal and 
double-buoyancy convection is the appearance of an additional 
"double-diffusive" eigenmode. The remaining two modes are 
analogous to the two possible modes of purely thermal con- 
vection, as figure [3] clearly demonstrates. In figure [3] and in 
the following, these three possible modes are denoted by aDD 
(additional Double-Diffusive mode), uTL (unstable Thermal- 
Like mode), and sTL (stable Thermal-Like mode). The aDD 
mode becomes unstable for smaller values of R, compared to 
the uTL mode. This provides a possibil i ty for l ow-Rayleigh 
number convection as suggested bv lBussd (12002bl) . The growth 
rate of the aDD mode is a non monotonic function of R,, and it 
is remarkable that in the case of a destabilizing compositional 
gradient (R s > 0), the aDD mode regains stability before the 
uTL mode becomes unstable. This limits the parameter space 
where low-Rayleigh number convection occurs, and indicates 
the existence of isolated regions of instability. 

The regions of linear stability (cr < 0) and instability (cr > 0) 
in the parameter space are separated from each other by a neu- 
tral surface. It is defined in implicit form by the characteristic 
equation of the eigenvalue problem © where cr = 0, 



(iu)P + a 2 )(i(DP + a 2 /L)((iu) + a 2 )a 2 + iarj*) 

-a 2 R,(i(A)P + a 2 /L) - a 2 R s (ia>P + a 2 ) = 0. 



(9) 



Following iBussel ( l2002bl) . this equation is split into real and 
imaginary parts from which the frequency of oscillations a> and 
the critical value of any parameter of the problem as a func- 
tion of the remaining ones can be found in explicit analytical 
form. Here, the thermal Rayleigh number R, is chosen as the 
principal control parameter, because it offers the possibility of 
direct comparison with the well-studied purely thermal case. 
The five-dimensional neutral surface, R, = R t (P, rf, L,R S , a), is 
represented graphically by its projections (neutral curves) onto 
the planes a - R,, R s - R,, rf - R,, and P - R, in panels (a, 
b, c, d) of figure ^respectively, for fixed values of the remain- 
ing parameters of the problem. The dependence on the Lewis 
number L is shown in the form of contour lines in this figure 
thus exhausting all possible parameter dependencies. The same 
approach is adopted to represent the corresponding frequency 
co = co(P,rf ,L,R s ,a) in figure [5] where |w| is plotted instead 
to allow visualisation of finer details in the plots. The most 
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Figure 3: (Color online.) The growth rates, <x = Re(/1), of the eigenmodes of double-diffusive convection in the rotating annulus geometry as a function of the 
thermal Rayleigh number R, at P = 1, tj* = 400, a = 4, L = 20, R s = 500 (solid red line) of R s = -500 (dashed blue line). The three possible modes are denoted 
by aDD, uTL, and sTL in the example of R s = 500. The basic state of no flow is linearly stable in the shaded region where the growth rate i s negative. Convection 
sets in the non-shaded region. The green dotted lines correspond to the well-studied purely thermal modes of convection e.g. (Busse, 1986) at the same parameter 
values (and R s = 0, L - arbitrary). The right panel is identical to the left one, only the scale of the y axis is enlarged to show finer details. 



prominent feature of the neutral curves is that they are multi- 
valued, and may split into closed, entirely isolated branches. 
This can be understood from the fact that the dispersion rela- 
tion ([9jl is a linear equation in R t , and a cubic equation in a>, so 
it has either one, two or three real roots as its discriminant takes 
negative, zero and positive values when parameter values are 
continuously varied. The stability of the basic state in the vari- 
ous regions formed thereby can be determined from the sign of 
the growth rate cr as described in relation to figure[3] For exam- 
ple in the case L = 30 of figure |4] convection occurs within the 
regions that have been shaded. 

The topology of the neutral curves of double-diffusive con- 
vection is essentially different from that in the case of purely 
thermal convection, also shown in figure |4] While in the latter 
case for fixed values of the other parameters there is one and 
only one critical value of R, above which convection occurs, 
in the former case up to three values of R, are needed to spec- 
ify stability criteria due to the multi-valued nature of the neu- 
tral curves. Note that the critical wave numbers associated with 
each of the three distinct critical values of R, are also different as 
seen in figure[4f a). Occurrence of isolated regions of secondary 
instability has been reported i n the case of qu a si-geo strophic 
purely thermal convection by iPlaut and Bussd (120021) . Note, 
that this is quite different from the isolated regions of primary 
instability discussed in this paper. It is likely that the double- 
diffusive case will exhibit even more complex behaviour in its 
transition to tertiary states, and this will be subject for future 
study. Neutral curves with similar complex topology have been 
previously reported i n unrelated s ituations, e.g. a differentially 
heated inclined box ( Hard 1971), quiescent layers with den- 
sity dependent on t wo or more strat ifying agencies with dif- 
ferent diffusivities jPearlsteinl 1981 ), isothermal shear flows 
dMeseguer and Marquesi 12002), buoyanc y-driven flows in an 
inclined layer dChen and Pearlstein , 1989b. and penetra tive con- 
vection in porous media dStraughan and Walkeru 19971) . 



It is of interest to discuss the expressions 



1 / rfP 



.>(!) 



a 2 R 2 
rf 2 P 

aRs 

fl 2 (l +P) *" TJ*P' 



2 \l+P 

rfa 



2PR S 

1 Tp' 



(10) 



and 



n(2) _ i 

a 2 rf »2p 



a 6 a 2 



(id 



^ = - 



aRs 
~rf~P 



1 + 



a 2 R s (l+P) 



T]* 2 P 



derived by [Busse d2002bl) as solutions to the dispersion relation 
((9]) in the asymptotic limit of large L. The first ro ot correspond s 
to the well-studied thermal Rossby waves, (e.g. iBusse , Il986h . 
modified by the presence of the second buoyancy component 
and describes the onset of the uTL mode. The physical na- 
ture of the second root ("the slow mode") can be understood 
from the observation that in the limit of large 77* the second 
term in (fTTT > vanishes and the critical Rayleigh number for the 
onset of Rayleigh- Benard convec tion in a non-rotating plane 
layer is recovered dBussel l2002bb . Thus, the additional buoy- 
ancy provided by the compositional gradient, R s d y r, counter- 
acts the unbalanced part of the convection-inhibiting Coriolis 
force, rfdyij/, in equations (01. Expressions ( TTOb and (flTT i are 
shown in figures [4] and [5] and it can be seen that they provide a 
good approximation to some pieces of the neutral curves even 
for moderate an d small values o f L and if . It has been implic- 
itly assumed by Busse d2002bl) that there is a unique critical 
Rayleigh number above which convection sets in. This led to 
the conclusion that the slow mode is the one preferred at onset, 



as Rf 1 is always smaller than ^ iJ . The presented results show 
that this assumption is not always correct, and that the multi- 
valued nature of the neutral curves must be taken into account. 
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Figure 4: (Color online.) Neutral curves of double-diffusive convection in the rotating annulus geometry. Projections of the neutral surfaces onto (a) the a - R t 
plane, (b) the R s — R t plane, (c) the rf - R, plane, and (d) the P — R t plane. In all panels, the values of a = 5, P = 10, rf = 600, R s = -500 (thick dashed blue 
lines) and R s = 500 (thin solid red lines), and L = 17 (innermost contour), 20, 30, 40 are kept fixed, except where they are given on the abscissa. As an example, the 
linearly unstable regions are shaded in the case R s = 500, L = 30; the other curves form si milar regions as well. The thick dotted green lines (a single point in panel 
(b)) correspond to the well-known purely thermal Rossby wave modes of convection e.g. (Busse, 1986) at the same parameter values (and R s = 0, L- arbitrary), and 
approximate closely the first asymptotic root dlOI . The black dash-dotted line represents the second asymptotic root 1 11 11 . 




For example, when the concentration gradient is destabilizing, 
R t is, actually, the value at which convection decays as R, is 
increased. 

4. Double-diffusive convection at finite amplitudes 

The linear results of section [3]demonstrate that low-Rayleigh 
number convection is indeed possible albeit the situation is 
more complicated. Below, the question whether such low- 
Rayleigh number flows are sufficiently vigorous to generate 
magnetic field is addressed and finite-amplitude properties of 



double-diffusive convection are explored. Finite-amplitude so- 
lutions are characterized by their mean zonal flow, stream func- 
tion, temperature and concentration perturbations, defined as 

v (x, t) = (8 x ifr) = 3^o, ® (x, t) = <©>, T (x, t) = <F>, 

where (f(y)) = L~ l J Q f(y)dy and L v = 2n/fi is the basic 
periodicity length, and by the amplitude of convection 

N,,N X 

a2 = Z 

/=!,«=! 
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Figure 6: (Color online.) (a) The time-averaged amplitude of convection (A), as a function of the thermal Rayleigh number R, in the case P = 10, L = 20, rf = 600, 
and R s = 500 (red) and R s = -500 (blue). In the case R s = 500, the region where convection occurs is shaded. For values of R, outside this region A is a decaying 
function of R h and the decay has been followed to values of A smaller than 10~ 20 in all cases. The numbers shown near the onset of convection indicate the 
preferred wavenumber a in each case, (b) The amplitude of convection for P = 10, L = 20, rf = 600 and R, = 14050, R s = (Case I, red dash-dotted line), R, = 0, 
R s = 17000 (Case II, blue dashed line) and R, = R s = 9100 (Case III, green solid line). 



Figure |6j a) shows the time-averaged flow amplitude, (A),, of a 
sequence of cases with increasing value of the thermal Rayleigh 
number R, and fixed values of the remaining parameters. In 
full agreement with the linear theory, two regions of convection 
are found, labeled I\ and h in this figure. They are separated 
by a region of vanishing flow. The amplitude of convection 
in region l\ is more than an order of magnitude smaller then 
that of the flow in region h- Comparison with figure [3] indi- 
cates that the low-amplitude flow in I\ is associated with the 
aDD modes which are characterised by relatively small values 
of cr, while the high-amplitude convection in h is likely asso- 
ciated with the uTL modes. Because of its small amplitude, 
low-Rayleigh number double-diffusive convection in region I\ 
is unlikely to be able to generate and sustain magnetic fields 
on its own as will be further discussed below. Within region 
I\ all computed solutions are stationary, and for this reason not 
illustrated, while as R, is increased in region h a sequence of 
stationary, time-periodic, quasi-periodic and chaotic solutions 
similar to those described in previous studie s of purely thermal 
convection, e.g. dBrummell and Hartlll993l) . is observed. 

The additional physics introduced by the second buoy- 
ancy force makes it difficult to compare directly double- 
diffusive convection to the much-better studied purely ther- 
mal case. A meaningful approach for comparison is to con- 
sider cases with equally large amplitudes. This is suggested 
by self-consistent MHD dynamo simulations where it has 
been established that sufficiently vigorous turbulent flow is the 
primary c ondition for generation of self-sustained magnetic 
fields e.g. ( Simitev and Busse . 2005 : Kutzner and Christensenl 
2002). For a comprehensive comparison the amplitude of the 
flow as a function, for instance, of the thermal and composi- 
tional Rayleigh numbers need to be computed. Then a contour 
plot of the data A(R S ,R,) can be a useful comparison map as 
cases located on the same energy level are expected to have sim- 



ilar ability for magnetic field generation. However, the practical 
computation of such a surface has proven too expensive even 
for the relatively simple annulus model considered here. For 
this reason, the attention is restricted below to a comparison of 
three representative cases: a purely thermal case, a purely com- 
positional case, and a mixed double-diffusive case, henceforth 
Cases I, II and III, respectively. The time-averaged amplitudes 
of convection in Cases I, II and III are (A), = 0.42, 0.42, 0.43, 
respectively. Although the values are not strictly equal, addi- 
tional simulations suggest that such small differences in ampli- 
tude are not essential for the intended comparison. The three 
cases have destabilizing thermal and compositional gradients, 
which is thought to be appropriate for the Earth's core. Pur- 
posefully, the cases are moderately rather than strongly driven 
to illustrate how simple known properties are affected by the 
presence of a second buoyancy. At these amplitudes the flows 
considered are associated with the uTL modes discussed previ- 
ously, rather than with the newly-found aDD mode. This choice 
is justified as the aDD modes do not produce sufficiently vig- 
orous flows with interesting structure, as already discussed in 
relation to figure|6la). 

Figure|6jb) demonstrates that for comparable time-averaged 
amplitude, the purely compositional Case II and the mixed Case 
III have a highly chaotic time dependence while the purely ther- 
mal Case I is stationary. The spatial properties of convection 
are shown in figure Q where the streamlines of the flow are plot- 
ted for the three cases along with the fluctuating parts of the 
temperature perturbation - ©o and the concentration pertur- 
bation F - Fo. The plots represent snapshots at fixed moments 
in time but they have been found to be representative for the 
three cases. The purely thermal Case I shows a regular roll- 
like pattern which does not change in time, while the struc- 
tures corresponding to Cases II and III are irregular and no pe- 
riodic behaviour of the patterns in time can be detected. These 
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Figure 7: (Color online.) The non-axisymmetric parts of the streamlines ip -"Vq = const, (contour lines, first row), the temperature perturbation - ©o (density 
plot, second row) and the concentration perturbation r - Tq (density plot, third row). The first, second and third columns correspond to the purely thermal Case I, 
to the purely compositional Case II, and to the double diffusive Case III, described in the caption of figure|6}b). 



differences may be explained by the fact that Case I with ra- 
tio Ri/Ri.dit = 1.52 is far less supercritical than Case II where 
RJRs.cnt = 17.84. The predominant wave number of convec- 
tion appears to be the same in all three Cases, and remains equal 
to 7 throughout the simulations. In comparison with the temper- 
ature perturbation that shows relatively broad roll structures, the 
concentration perturbation forms thinner plume-like structures, 
consistent with the smaller compositional diffusivity. The time- 
and azimuthaly-averaged properties of convection in the three 
cases are compared in figure [8] The most obvious difference is 
observed in the profiles of the time-averaged mean zonal flow 
and the Reynolds stress. These quantities are, indeed, related in 
that t he mean flow is generated primarily by the Reynolds stress 
(e.g. IPlaut et all 120081: iBussei l2002bl) . While in the purely 



thermal Case I the mean flow is symmetric with respect to the 
mid-channel x — 0, and retrograde at its ends x = +1/2, in the 
purely compositional Case II it is asymmetric with respect to 
x = 0, retrograde at x = 1/2 and prograde at x = -1/2. This 
asymmetry can be explained by the property that, unlike in the 
purely thermal case, the value of R s in the com positional case 
is bey ond the onset of the mean-flow instability (lOr and Busse , 



1987). The mean flow in the mixed Case III appears similar to 



the purely compositional case. The remaining panels in figure 
[8] show that the mean properties of the mixed case are simi- 
lar to the corresponding ones of the pure cases. In summary, 
it appears that double diffusive convection associated with the 
uTL modes can be understood on the basis of the correspond- 
ing single-diffusive cases, and that purely-thermal convection 
is more efficient in imprinting its properties on the overall flow 
even when less supercritic al. This conclusion is co nfirmed by 
the experimental results of ICardin and Olsonl (119921 ) who stud- 
ied thermochemical convection in rotating spherical shells and 
found that the structure of thermochemical flows is more like 
that of purely thermal convection. 



5. Conclusion 

Convection driven by density variations due to differences 
in temperature and concentration diffusing at different rates in 
a rotating cylindrical annulus with conical end caps has been 
studied. It is shown by a linear analysis that the neutral surface 
describing the onset of convection in this case has an essentially 
different topology from that of the well-studied purely ther- 
mal case. In particular, due to an additional "double-diffusive" 
eigenmode (aDD), neutral curves are typically multi-valued 
and form regions of instability in the parameter space which 
may be entirely disconnected from each other. It is confirmed 
that the asymptotic expres sions for the cri tical Rayleigh number 
and frequency derived by IBussei (12002bl) describe the onset of 
convection over an extended range of non-asymptotic param- 
eter values but do not capture the full complexity of the neu- 
tral curves. The results necessitate a revision of the assumption 
that there is a unique critical value of the control parameter, 
e.g. R,, and call for a better appreciation of the multivalued na- 
ture of the critical curves. It is been found that finite-amplitude 
low-Rayleigh number convection due to aDD modes is possible 
over a wide parameter range. However, the resulting flow am- 
plitudes are significantly lower than those of due to the familiar 
uTL modes of convection. For this reason, low-Rayleigh num- 
ber flows are unlikely to be able to generate and sustain mag- 
netic fields on their own. In order to address a more geophys- 
ically relevant situation, the nonlinear properties of convec- 
tion are then investigated in the case when both driving agen- 
cies are destabilizing and produce sufficiently vigorous flow. It 
is proposed that a meaningful approach for direct comparison 
of finite-amplitude double-diffusive convection and the better 
studied single-diffusive case is to compare flows with equally 
large kinetic energies. Using this criterion the characteristics of 
a purely thermal case, a purely compositional case and a mixed 
driving case are compared. As similar flow amplitudes purely 
compositional and double-diffusive cases are characterized by 
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<i'o>, (vu)t <©o>, (To), 



Figure 8: (Color online.) Profiles of the time-averaged (a) mean velocity (vg)t, (b) Reynolds stress (vu) t , (c) mean temperature perturbation (&o)t, and (d) the mean 
concentration perturbation <ro> f . Red dash-dotted lines indicate Case I, blue dashed lines indicate Case II, and green solid lines indicate Case III, described in the 
caption of figure|6jb). 



a stronger time dependence compared to purely thermal cases, 
and by prograde mean zonal flow near the inner cylindrical sur- 
face. It is argued that double-diffusive cases may be understood 
on the basis of purely driven ones. 

Although, its low amplitude is likely to prevent double- 
diffusive convection at values of the Rayleigh number signif- 
icantly lower than those for single-diffusive convection from 
generating magnetic fields in the bulk of planetary cores, it 
is tempting to speculate that this type of flow may have im- 
portant effects in stratified layers located just under the core- 
mantle boundary. Several mechanisms have been suggested 
for the for the possible formation such layers, including the 
build-up of light elements released during inner core solidi- 
fication (iBraginskyl 120061) . and thermal or chemical interac- 



tion between the mantle and the core ( Fearn and Looer, 11981, 
Lister and Buffeti 1 19951 iBuffett and SeagleL l2010h . Crucially. 



evidence for stratification has been recently reported in seis- 
mic o bservations of the outer core (IHelffrich and Kaneshima , 
20101) . Models of inert stably stratified outer layers have been 
found to produce magnetic fields with morphology rather dis- 
similar to that of the observed field because of a thermal 
wind that produces unfavorab l e zon al flows throughout the 
core ([Stanley and Mohammadii 120081) . Inert layers, have also 
been found to behave like a no -slip virtual boundary for the 



convective motion underneath ( Takehiro et al. , 2010b . This 



last finding imposes a significant constraint on the flow, as 
it is well known that convection structures and the morphol- 
ogy of the magnetic field crucially depend on the boundary 
conditions ( Simitev and Busse, 2005tlKutzner and Christensen , 



2002 ; Sakuraba and Roberts , 2009b . The situation may be sig 
nificantly different if the stratified layer is convecting (even 
weakly) rather than inert and the low-Rayleigh number regime 
I\ found here offers one such possibility. This possibility will 
be subject of future research. In addition, it will be of interest to 



investigate whether the results reported in this paper hold in the 
more realistic case of a spherical shell. In particular, the spher- 
ical case may allow the aDD modes to grow to a much larger 
amplitude, because geostrophy is not hard-wired into the for- 
mulation of the spherical model as it is in the annulus case. If 
this should be the case, low Rayleigh-number convection may 
have a more significant role in core dynamics. The influence 
of imposed magnetic fields and the general parameter depen- 
dences of the problem must also be studied in more detail to 
explore scaling relationships and the possibility of further in- 
teresting dynamics. 
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